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> : ABSTRACT 
(N 

. Within the last 10 years, strong, narrow Na I absorption has appeared at Vq ~ —38 

', km s~^ toward the halo star HD 219188; that absorption has continued to strengthen, 
by a factor 2-3, over the past three years. The line of sight appears to be moving 

fH _ into/through a relatively cold, quiescent intermediate velocity (IV) cloud, due to the 

Qh' 13 mas/yr proper motion of HD 219188; the variations in Na I probe length scales of 

O ■ 2-38 AU/yr. UV spectra obtained with the HST GHRS in 1994-1995 suggest iV(Htot) 

^ ■ ~ 4.8 X 10^"^ cm~^, "halo cloud" depletions, nn ~ 25 cm~^, and rig ~ 0.85-6.2 cm~^ (if 

Ci . r ~ 100 K) for the portion of the IV cloud sampled at that time. The relatively high 
fractional ionization, Ue/nu > 0.034, implies that hydrogen must be partially ionized. 

^ I The A^(Na l)/A^(Htot) ratio is very high; in this case, the variations in Na I do not imply 

H ' large local pressures or densities. 

Subject headings: ISM: clouds — ISM: structure — Line: profiles — stars: individual 
(HD 219188) 
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1. Introduction 

Interstellar absorption lines were first recognized as interstellar by virtue of their invariance 
amidst the shifting stellar lines in the S Ori A spectroscopic binary system (Hartmann 1904). Re- 
cently, however, it has become evident that some clearly interstellar lines do vary in strength and/or 
velocity. Variations in Na I and/or Ca II have been observed toward several stars in the direction 
of the Vela supernova remnant (Hobbs ct al. 1991; Danks & Scmbacli 1995; Cha &: Sembach 
2000) — presumably due to cloud motions in that energetic environment. High A^(Ca Il)/A^(Na I) 
ratios suggest significant shock processing of the interstellar dust grains in some of those clouds. 
In several other cases, however, the narrow Na I line widths (6 ~ 0.3-0.5 km s"-*^) and/or low 
iV(Ca II)/iV(Na I) ratios characterizing the variable components suggest that the gas is relatively 
cool and quiescent (Blades et al. 1997; Price et al. 2000; Crawford ct al. 2000; Lauroesch, Meyer, &i 
Blades 2000). Detailed information regarding the physical conditions in these varying components 
is generally not available, however. 

These indications of temporal variability may be related to recent results on small scale spatial 
structure in the ISM. Observations of H I 21cm absorption toward high-velocity pulsars have 
revealed structure on scales of 5-100 AU (Frail et al. 1994). Observations of absorption due to the 
trace ions Na I and/or K I toward globular clusters have found significant differences on scales of 
10^-10^ AU (Langer, Prosser, 8z Sneden 1990; Kemp, Bates, & Lyons 1993; Meyer & Lauroesch 
1999). Differences in Na I and/or K I absorption toward the individual members of binary or 
multiple stellar systems are commonplace, suggesting ubiquitous structure on intermediate scales 
of lO^-lO'' AU (Watson & Meyer 1996; Lauroesch & Meyer 1999). If these differences are due to 
variations in overall hydrogen column density [assuming "typical" A'^(Na I)/A^(H) ratios], the high 
implied densities (generally several thousand cm~^) are difficult to reconcile with clouds in thermal 
pressure equilibrium at typical interstellar pressures of order 3000 cm~^K. The differences may 
reflect a population of cold, dense filaments or sheets, embedded in warmer, less dense neutral gas 
and containing 10-30% of the total column density of cold, neutral gas (Heiles 1997). Alternatively, 
the differences in the trace ions might be due to variations in density and/or ionization on those 
small scales (e.g., Lauroesch et al. 1998). Differences in ionization would not explain the variations 
seen in H I, however. 

HD 219188 is an early B supergiant located about 2800 pc from the sun at {l,b) ~ (83°,— 50°). 
Because of its location and brightness {V ~ 6.9), it has been observed in a number of studies 
investigating the interstellar clouds in the lower Galactic halo (Miinch & Zirin 1961; Albert 1983; 
Danly 1989; Sembach, Danks, & Savage 1993). Strong, multi-component absorption from Na I and 
Ca II is apparent for —17 km s^^ ^ ^0 ^ +5 km s""*^; weaker Ca II absorption extends to about 
—35 km s~^. High resolution Na I spectra obtained in 1997, however, revealed a new, relatively 
strong, narrow component at f ~ —38 km s~^. This component is also present in Ca II and in 
several species in UV (GHRS echelle) spectra obtained in 1994-1995. Subsequent monitoring of the 
line of sight shows that the column densities of both Na I and Ca II have increased by additional 
factors of 2-3 since 1997. In the following sections, we discuss the recent optical and UV spectra 
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of HD 219188 and the inferred properties of the newly-revealed intermediate-velocity cloud along 
that line of sight. 

2. Optical Data 

The initial Na I and Ca II spectra of HD 219188 were obtained with the coude feed and coude 
spectrograph at the Kitt Peak National Observatory, as part of a program to acquire high resolution 
(FWHM < 1.5 km s~^) optical spectra of halo stars previously observed at slightly lower resolution 
in the UV with the HST GHRS (see Fitzpatrick & Spitzer 1997). After the 1997 Na I spectrum 
of HD 219188 revealed the new intermediate- velocity (IV) component at —38 km s"^, additional 
spectra were obtained during various runs at Kitt Peak, the Anglo-Australian Observatory, and 
the European Southern Observatory in order to monitor possible further variations (Table 1). 

The spectra were processed and extracted using standard procedures within IRAF, then nor- 
malized via polynomial fits to continuum regions adjacent to the absorption lines. The resulting 
line profiles were fitted with multiple Voigt profile components, yielding column densities, line 
widths, and heliocentric velocities for the individual components (e.g., Welty, Morton, & Hobbs 
1996). Some of the Na I and Ca II profiles are shown in Figure 1, with the earlier profiles described 
by Albert (1983) and Sembach et al. (1993) at the top for comparison. 

The changes in the strength of the —38 km component — over the past 20 years and 
over intervals as short as 6 months — are quite evident and striking (Fig. 1; Table 1). While 
weak Na I absorption at u ^ —38 km appears to have been present in 1991 (Sembach et al. 

1993), none was detected in 1980 (Albert 1983). By 1997, the Na I column density was roughly a 
factor of 10 higher than the limit in Albert (1983), and it has increased by a further factor of 2 
since then. As corresponding IV Ca II absorption may be present in the profiles from both earlier 
studies (Fig. 1), we have fitted those profiles using the component velocity structure derived from 
our higher resolution spectra to estimate A/"(Ca II) in the IV cloud at those earlier epochs. While 
our Ca II column densities in 1995 and 1997 may be consistent with the (uncertain) value reported 
by Albert (1983), A^(Ca II) has increased by a factor of 2-3 since then. Over the past 3-5 years, 
the A^(Na l)/A^(Ca II) ratio thus appears to have remained roughly constant (though significantly 
higher than previous values or limits) (Table 1). 

In contrast to the changes in column density, neither the velocity of the IV component nor the 
Na I line width has changed significantly since 1995. The width of the Na I line, b ~ 0.55-0.60 km 

implies a maximum temperature for the cloud of 490 K and a maximum 3-dimensional internal 
turbulent velocity (3^/^)wt = 0.73 km s~^ (comparable to the sound speed for T ~ 80-100 K). 
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3. UV Data 

HD 219188 was observed with the HST GHRS on 1994 June 6 and 1995 May 16-17, at 
resolutions of about 3.5 km s~^, under GTO program 5093. Twelve wavelength settings were 
obtained, using the FP-SPLIT procedure to enable the identification and ehmination of detector 

fixed-pattern features. The individual sub-exposures for each wavelength setting were combined 
using the POFFSETS and SPECALIGN routines within the STSDAS package. The summed 
spectra were normalized via polynomial fits to continuum regions adjacent to the absorption lines; 
the empirical S/Ns range from 30-50. The interstellar absorption line profiles were then fitted using 
the component structures obtained from the higher resolution spectra of Na I and Ca II (Fitzpatrick 
& Spitzer 1997; Welty et al. 1999). 

The —38 km component is detected in the ground state lines of C I, C II, O I, Si II, S II, 
and Fe II and in the excited fine-structure lines of C I*, C I**, and C II* (Table 2, Fig. 2). If 

sulfur is undepleted, then A^(S II) suggests a total IV hydrogen column density A''(Htot) = A^(H I) 
-|- A^(H II) ~ 4.8 X 10^''' cm~^. The IV column density estimated for O I (which should track H I) 
is consistent with that A^(Htot) and the typical interstellar oxygen abundance [0/H] ~ —0.4 dex 
(Meyer, Jura, & Cardelli 1998) — which suggests that the IV gas is largely neutral. The C II A1334 
line is saturated, but we may use A^(S II) to estimate the total IV carbon abundance [N{C 1) -|- 

A^(C II)] ~ 6.8±2.1 X 10^3 cm-2, as [C/S] 0.4±0.1 dex in diffuse Galactic clouds. The column 

densities for Fe II and Si II yield depletions similar to those found for clouds in the Galactic halo 
(Welty et al. 1999, and references therein). 

Analysis of the C I fine-structure excitation equilibrium enables estimates for the thermal 
pressure and local density in the IV gas (Jenkins & Shaya 1979). The observed A''(C I*)/A^(C I) = 
0.17±0.03 and N{C I**)/iV(C I) = 0.06±0.02 yield n^T ~ 2500 cm^^K and nn ~ 25 cm-^ for T 
= 100 K. The thickness of the region probed, A'^(H)/nH, would thus be ~ 0.006 pc, or about 1300 
AU. Because n^T and nn depend on the assumed T, we have listed values for several temperatures 
consistent with the limit from 6(Na I) in Table 3. 

Estimates for the electron density in the IV gas may be derived both from considerations of 
photoionization equilibrium and from C II excitation equilibrium. If photoionization equilibrium 
is assumed for carbon [with V/a = 24 for the WJl interstellar radiation field and T = 100 K 
(Pequignot & Aldrovandi 1986) and with N{C II) as estimated above], then N{C l)/N{C II) = Ue 
(r/a)~^ yields ne(C I) ~ 6.2ib2.4 cm~^. Analysis of the C II fine-structure excitation (Fitzpatrick 
& Spitzer 1997; Welty et al. 1999) for T = 100 K, however, yields ne(C II) ~ 0.85+[j:|^ cm'^ — 
about a factor of 7 smaller. While the uncertainties in both Ue estimates are dominated by the 
almost 50% uncertainty in iV(C II), the ratio of the two estimates is not very sensitive to N{C II). 

Similar discrepancies between ne(C I) and ne(C II) have been noted for other lines of sight 
(Fitzpatrick k, Spitzer 1997; Welty et al. 1999). While C I might be concentrated in a smaller, 
denser region than C II, ne(C I) for that smaller region would be even larger than the average 
value derived above, and would imply a fractional ionization ne/nn greater than 25%. The ne(C I) 
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would be reduced, however, if the radiation field were weaker or if T were lower than assumed 
above. Alternatively, the large rig (relative to the assumed WJl radiation field) suggests that 
enough negatively charged large molecules (e.g., PAH's) could be present for charge exchange 
reactions (X+ + LM~ — > X*^ + LM*^) to significantly enhance the column densities of the trace 
neutral species — in which case ne(C I) would be too high (Lepp & Dalgarno 1988; Lepp et al. 
1988; Welty & Hobbs 2001). Even if that is the case, however, ne(C II) still yields ne/nn ~ 0.034 
— more than a factor 240 greater than the value 1.4 x lO^'^ that would be due to photoionization 
of carbon (usually assumed to be the primary source of electrons). Hydrogen therefore appears to 
be partially (> 3.4%) ionized in the —38 km s~^ cloud. It is not clear how such a high fractional 
ionization could be maintained. 



4. Discussion 

The Leiden-Dwingeloo H I maps (Hartmann & Burton 1997) suggest that the line of sight to 
HD 219188 lies near the inner edge of a shell of IV gas, seen over the range —30 to —60 km s~^ 
(^0). Despite its possible association with an expanding and/or compressed shell, the cloud's fairly 
low T, fairly high iV(Na l)/A^(Ca II) ratio, moderate pressure, and stable velocity suggest that it is 
relatively quiescent. The variations in absorption may instead be due to the relatively high proper 
motion of HD 219188 about 13.4 mas/year (ESA 1997) — which corresponds to a transverse 
velocity of about 180 km s^^, or about 38 AU per year, at 2800 pc. The steadily increasing column 
densities of Na I and Ca II and the high fractional ionization suggest that the line of sight is moving 
into/through the partially ionized edge of the IV cloud due to the motion of the background star. 

Unfortunately, the distance to this IV cloud is not well constrained. A lower limit of 150 pc 
is provided by the lack of IV absorption toward (3 Psc (<2°away, but toward detectable IV H I 
emission); the upper limit (2800 pc) is set by HD 219188 itself. The length scale in the IV cloud 
probed by the stellar proper motion is thus between 2 and 38 AU/ycar — considerably smaller 
than the 10^-10^ AU scales sampled by the binary star studies. The cloud distance (and thus the 
length scale) could be better constrained by obtaining Na I spectra of other stars in this region (at 
various distances). 

While we do not have a spectrum of Na I from 1994-1995, when the UV spectra were obtained, 
estimates based on the C I and Ca II observed then suggest that A'^(Na I) could have been ~ 2 
X lO^-*^ cm~^. The ratio iV(Na I)/iV(H) ~ 4.2 x 10"^ would thus have been much higher than is 
typically observed for diffuse clouds with iV(H) > 10^^ cm~^ (Welty & Hobbs 2001) — prcsTmiably 
due to the high fractional ionization. Furthermore, if Na is depleted by a factor 4 in the IV gas 
(Welty k, Hobbs 2001), then Na I would have been the dominant form of Na. 

If such high iV(Na I)/iV(H) ratios and fractional ionizations are representative of other sight- 
lines where differences in Na I have been observed over small transverse scales, then the corre- 
sponding differences in iV(H) — and the resulting inferred pressures and densities — would not 
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be as large as previously estimated. Lauroesch et al. (2000) have estimated ne/nu ^ 0.008 for 
the variable component at u ~ 9 km toward HD 32040, and N(Na l)/A^(H) is likely high [for 
A^(H) ~ 6 X 10^^ cm~^] in the —8.6 km s"^ component toward fi^ Cru (Meyer & Blades 1996). 
The high A'^(Na I)/A^(H) ratios found for a sample of high- and intermediate- velocity clouds with 
A'^(H) < 10^^ cm~^ (Wakker & Mathis 2000) may also generally be due to partial ionization of the 
hydrogen. These lower N(}1) clouds may thus represent a different regime from the more familiar 
diffuse neutral clouds with Ar(H) > 10^^ cm^^. 

Lauroesch et al. (1998) have conjectured that some of the small-scale structure observed 
in Na I may be due to fluctuations in ionization [rather than in overall Ar(H)] for cases where 
corresponding differences are not seen for Zn II [which should track Ar(H)]. The low overall column 
densities, moderate local densities, and high fractional ionization found for the IV cloud toward HD 
219188 suggest that some of the small-scale structure seen in Na I and K I may be due to partially 
ionized, relatively low-density cloud edges and/or fragments — rather than dense, neutral clumps. 
Further monitoring of the HD 219188 line of sight over the next several years — particularly in 
the UV with HST/STIS — should provide interesting information on the variations of density, 
temperature, and depletions with depth in the IV cloud. 

We thank D. Willmarth (KPNO), S. Ryan (AAO), and M. Kucrstcr (ESO) for their help in 
setting up the spectrographs, Lew Hobbs for pointing out the large proper motion of HD 219188, 
and Bart Wakker for providing a map of the 21 cm emission (and individual spectra) near HD 
219188. DEW acknowledges support from NASA LTSA grant NAGW-4445. 
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Table 1. Optical Data: v = —38 km s Component^ 



Date 


Observatory 
(Instrument) 


Resolution 
(km s-i) 


(cm-2) 


Na I - - 

b 

(km s-i) 


V 

(km s-i) 


(cm-2) 


- - - Ca II - 
b 

(km s-i) 


V 

(km s-l) 


Ar(Na I)/Ar(Ca II) 


1980.75"= 


McDonald (2.7m) 


5.9/5.4"= 


<3.0 






4±1 


(1.5) 


(-38.6) 


<0.8 


1991.56<^ 


ESO (CAT/CES) 


4.4 


4±3 


1.5±0.9 


-37.9±1.1 


2±1 


(1.5) 


(-38.6) 


2.0±1.8 


1995.82 


KPNO (coude feed) 


1.40 








3.1±0.3 


(1.5) 


-38.6±0.2 




1997.77 


KPNO (coudc feed) 


1.35 


27±1 


0.59±0.06 


-38.3±0.1 


4.9±1.0 


(1.5) 


-38.9±0.3 


5.5±1.1 


1998.68 


AAO (AAT/UCLES) 


5.0 


34±2 


(0.6) 


-38.0±0.1 










1999.42 


KPNO (coude feed) 


1.50 


39±4 


(0.6) 


-38.1±0.1 










1999.98 


ESO (3.6m/CES) 


1.20 


52±2 


0.55±0.05 


-38.2±0.1 










2000.46 


KPNO (coude feed) 


1.50 


60±5 


0.56±0.09 


-38.3±0.1 


10.7±1.2 


(0.6) 


-38.3±0.1 


5.6±0.8 



^Uncertainties are 1<t; values in parentheses were fixed in the fits. 
''Column density in units of 10^*^ cm^^. 

"= Albert 1983; FWHM for Na I/Ca II; A'(Ca II) from new fit to generated spectrum. 

"^Sembach et al. 1993; shifted to align low-velocity absorption; N{Ca, II) from new fit to generated spectrum. 
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Table 2. UV Data: v = -38 km s'^ Component^ 



Ion 


Line 

(A) 


N{X) 
(cm-2) 


[X/S]^ 


C I 


1560 


10.7±1.0el2 




C I* 


1560 


2.3±0.3el2 




C I** 


1561 


0.8±0.3el2 




C II 


1334 


(4.9±2.1el3)^ 


(-0.4)^ 


C II* 


1335 


5±lel2 




I 


1302 


l.l±0.4el4 


-0.5 


Si II 


1304 


9±lel2 


-0.3 


S II 


1253 


9±2el2 




Fe II 


2374, 2600 


4.0±0.4el2 


-0.6 



^Uncertainties are la; limits are 2a. 

^[X/S] = log{A^(X II)/A^(S II)} - 
log{ (Ax/As)©}, using the solar system 
meteoritic abundances of Anders & Grevesse 
1989 and Grevesse & Noels 1993. 

'^Estimated assuming typical interstellar 
[C/S] = -0.4±0.1 dex (for C I + C II). 
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Table 3. Inferred Properties: v = —38 km s Component 



Quantity 




Value 




N{B.) (cm-2) 




4.8 X 10^^ 




r(K) 




<490 






25 K 


100 K 


400 K 


log(nHr) (cm-^K) 


3.4±0.1 


3.4±0.1 


3.8±0.1 


riH (cm~^) 


100 


25 


15 


thickness (AU) 


325 


1300 


2200 


ne{C I)^(cm-3) 


2.6 


6.2 


14.7 


ne(C II)'^(cm-3) 




0.85 


0.79 






0.034 


0.053 



^Assuming photoionization equilibrium for car- 
bon; WJl radiation field. 

''From C II fine-structure excitation equilibrium; 
no valid value for T = 25 K. 

^Using He from analysis of C II fine-structure. 
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Fig. 1. Selected Na I A5895 (left) and Ca II A3933 (right) spectra toward HD 219188. The sources 
and dates of the spectra are indicated; see Table 1. The column densities in the intermediate- velocity 
component at u ~ —38 km s~^ have increased by factors of 2-3 over the past 3-5 years; no Na I 
was detected at that velocity in 1980. The lower resolution spectra from Albert (1983; Alb80) and 
Sembach et al. (1993; SDS91) were generated from the component structures listed therein. 
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Fig. 2. — Profiles of selected UV lines toward HD 219188, observed with GHRS at resolutions of 
about 3.5 km in 1994-1995. The vertical dotted line indicates the component at —38 km s~^. 



